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Myostatin induces mitochondrial metabolic alteration 
and typical apoptosis in cancer cells 

Y Liu^'2'^ H Cheng^ ^ Y Zhou^'^ Y Zhu^ R Bian^ Y Chen^ C Li^ Q Ma^ ^ Q Zheng^ ^ Y Zhang^ '^ H Jin\ X Wang\ Q Chen*'^ "^ 
and D Zhu*'^'^ 

Myostatin, a member of the transforming growth factor-)? superfamily, regulates the glucose metabolism of muscle cells, while 
dysregulated myostatin activity is associated with a number of metabolic disorders, including muscle cachexia, obesity and type 
II diabetes. We observed that myostatin induced significant mitochondrial metabolic alterations and prolonged exposure of 
myostatin induced mitochondria-dependent apoptosis in cancer cells addicted to glycolysis. To address the underlying 
mechanism, we found that the protein levels of Hexokinase II (HKII) and voltage-dependent anion channel 1 (VDAC1), two key 
regulators of glucose metabolisms as well as metabolic stress-induced apoptosis, were negatively correlated. In particular, 
VDAC1 was dramatically upregulated in cells that are sensitive to myostatin treatment whereas HKII was downregulated and 
dissociated from mitochondria. Myostatin promoted the translocation of Bax from cytosol to mitochondria, and knockdown of 
VDAC1 inhibited myostatin-induced Bax translocation and apoptosis. These apoptotic changes can be partially rescued by 
repletion of ATP, or by ectopic expression of HKII, suggesting that perturbation of mitochondrial metabolism is causally linked 
with subsequent apoptosis. Our findings reveal novel function of myostatin in regulating mitochondrial metabolism and 
apoptosis in cancer cells. 
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Myostatin, a critical myokine and a member of the transform- 
ing growth factor-^^ (TGF-p) superfamily, acts as a negative 
regulator of muscle mass"*'^ and its mutation results in 
muscular hypertrophy and hyperplasia.^ "^ Like most Activin/ 
TGF-j5 family members, myostatin binds to the cell surface 
activin receptor II or MB (ACTRII, ACTRIIB), which recruits 
type I receptor ALK 4/5 (Activin receptor-like kinase 4 or 5) to 
form a complex.^'^ ALK5, which is also involved in the TGF-j5 
signaling pathway, could phosphorylate and activate Smad2/ 
3 for its nuclear translocation and regulation of target genes 
transcription.'' Myostatin also utilizes the non-canonical path- 
way, such as activation of the MARK pathway, or inhibition of 
the PI3K-Akt/GSK pathway, leading to suppression of 
myoblast proliferation and differentiation.^'^ Recent studies, 
including ours, have demonstrated that myostatin regulates 
glucose metabolism by promoting glucose consumption and 
uptake, increasing glycolysis and inhibiting glycogen synth- 
esis in skeletal muscle cells.''^''''' 

Myostatin circulates in the blood and its receptors are 
ubiquitously expressed in all tissues. Emerging evidence has 



suggested its function in regulating energy metabolism in both 
muscle and non-muscle cells. Knockout of myostatin in 
genetic mouse models of obesity and diabetes improved 
glucose metabolism and reduced obese phenotype."*^ More 
specifically, it was found that myostatin treatment inhibited 
glucose uptake in placental cells. ""^ Despite these tantalizing 
results, it is possible that the reduction in adipose tissue mass 
in myostatin mutant mice is an indirect result of metabolic 
changes in skeletal muscle. It remains to be explored 
whether and how myostatin regulates metabolism in non- 
skeletal muscle tissues. 

Accumulating evidence has also demonstrated that dysre- 
gulated myostatin is associated with metabolic disorders such 
as cachexia induced by tumors.''^ ''^ As most cancer cells 
express myostatin receptors and several members of Activin/ 
JGF-p family play very important roles in regulating cell 
growth, metabolism and apoptosis,^'' it is therefore concei- 
vable to hypothesize that myostatin exerts functional roles in 
regulating cancer cell growth or death by regulating energy 
metabolism. This is important since rapidly growing tumor 
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cells typically display altered aerobic glycolysis (Warburg 
effect)''^ and metabolic dysregulation is related to tumor 
growth and cell death. In this report, we therefore tested this 
hypothesis and our results demonstrated that myostatin 
induces metabolic shift from oxidative phosphorylation 
(OXPHOS) to glycolysis in cancer cells and interestingly the 
chronic exposure of myostatin results in the activation of 
mitochondria-dependent apoptosis. In an effort to understand 
the underlying mechanism, we showed that upregulation of 
VDAC1 and Bax translocation to the mitochondria played 
critical role in myostatin-induced apoptosis in cancer cells. 
The findings presented in this study suggest that a better 
understanding of the mechanistic links between cancerous 
metabolism and growth control by myostatin may be useful for 
developing better treatments of human cancer. 



Results 

Myostatin induces mitocliondrial metabolic alterations 
in cancer cells. To test the possibility that myostatin 
regulates mitochondrial metabolic activities in cancer 
cells, we first examined glucose consumption and lactate 
production in HeLa cells following treatment with myostatin. 
Glucose consumption (Figure 1a) and lactic acid production 
(Figure 1b) were significantly accelerated as early as 6 or 
12 h, respectively. A similar lactic acid production profile were 
also observed in various cancer cell lines including AZGY- 
83A cell (a lung adenocarcinoma cell line) and MGF-7 cell (a 
breast cancer cell line), but these changes were less drastic 
in SW480 (a colorectal cancer cell line) and SH-SY5Y (a 
neuroblastoma cell line) (Figure 1c), indicating the differential 
function of myostatin in regulating mitochondrial metabolism 
in cancer cells. In support of this observation, we also 
examined the metabolic fluxes of both glycolysis and 
OXPHOS, and not surprisingly, observed a significant 
increase of ECAR (extracellular acidic ratio; an indicator of 
glycolysis flux) whereas a robust decline of OCR (oxygen 
consumption ratio; an indicator of OXPHOS) from 12 to 24 h 
(Figures Id and e), which is supported by the observation 
that mstn-nuW mice showed slightly higher level of respiration 
ratio than their wild type littermates (Supplementary Figure 
SI). To further verify the effect of myostatin on mitochondrial 
metabolism, we assayed ATP concentration and found that 
the ATP level was maintained at early time points following 
treatment with myostatin, and that prolonged treatment 
caused drastic ATP reduction as early as 12h and almost 
complete depletion by 36 h (Figure If). In addition, the loss of 
mitochondrial membrane potential (A^m) became evident 
24 h after myostatin treatment as measured by DioCe,^ a 
commonly used dye to measure mitochondrial membrane 
potential (Figure 1g; Supplementary Figure S2). Collectively, 
these results demonstrate that myostatin alters mitochondrial 
energy metabolism by switching glucose metabolism from 
OXPHOS to glycolysis in glucose-addicted cancer cells. 

IVIyostatin induces apoptosis via mitochondrial apopto- 
tic pathway in HeLa cells. We observed that prolonged 
exposure of myostatin caused significant loss of viability of 
the cell. We thus tested if myostatin induced the loss of 



viability through a mechanism of apoptosis in HeLa cells. 
Annexin V-propidium iodate double staining clearly showed 
that myostatin induced phosphatidylserine externalization, a 
hallmark of apoptosis in a time- and dose-dependent manner 
(Figures 2a and b). Western blotting revealed a decrease in 
both procaspase-9 and procaspase-3 levels following treat- 
ment with myostatin for 24 h (Figure 2c). The cleavage 
product of PARP, a substrate of caspase-3, was evident after 
24 h of treatment (Figure 2c). We also measured the caspase 
activity in individual HeLa cells by FITC-VAD-FMK, a 
fluorescent marker specific for activated caspases. Gas- 
pases are activated after 24 h of myostatin treatment and can 
be potently inhibited by z-VAD, a pan-caspase inhibitor 
(Figure 2d). Furthermore, the presence of myostatin sig- 
nificantly reduced clonogenicity (Figure 2e), which further 
supports that myostatin induced apoptosis in HeLa cells. 
These results demonstrate that myostatin induced typical 
apoptosis, as characterized by the activation of caspase and 
the loss of clonogenicity. 

Myostatin induced drastic mitochondrial morphological 
changes from a filamentous shape into a donut shape 
(Supplementary Figure S3). Similar changes were observed 
in cells under metabolic stress-induced cell death. To assert 
that myostatin induces mitochondrial-dependent apoptosis, 
we assessed whether cytochrome c was redistributed from 
the mitochondria into the cytosol, where it activates caspase-9. 
Gytochrome c levels were significantly reduced in mitochon- 
drial fractions, while were increased in cytosolic fractions after 
24 h of treatment with myostatin, although the total cytochrome 
c level remained unchanged (Figures 2f and g). These results 
therefore showed that the mitochondrial pathway is involved in 
myostatin-induced apoptosis. 

Myostatin inhibits tumor growth in xenograft tumor 
models. We next examined whether myostatin induces 
tumor cell apoptosis in xenograft tumor models. We first test 
the effect of myostatin on xenograft tumors in immune- 
deficient nude mice. Treatment with varying concentrations 
of myostatin did not affect the body weight of the animals 
treated (data not shown). Notably, we observed substantial 
reduction of tumor weight in the mice treated with myostatin 
(Figures 3a and b). The average tumor size in the myostatin- 
treated group was significantly reduced, compared with that 
of the vehicle group (Figure 3c). As a matter of fact, after 
6 days of treatment with myostatin, 2 of the 10 treated mice 
achieved tumor regression below palpable levels. We then 
asked whether the reduced tumor size was due to tumor cells 
undergoing apoptosis as observed in in vitro cultured cancer 
cells. Indeed, TUNEL-staining analysis revealed a markedly 
increased number of apoptotic tumor cells in the myostatin- 
treated tumor compared with the apoptotic cells from vehicle 
(PBS)-treated tumors (Figures 3d and e; Supplementary 
Figure S4). Furthermore, functional role of myostatin in 
inhibiting tumor growth was also corroborated by approach- 
ing this question with an immune-competent xenograft 
mouse model. B16F10 cells, which are sensitive to myosta- 
tin-induced apoptosis (Supplementary Figure S5), were used 
to establish xenograft melanoma in G57BL/6 mice. Treat- 
ment with myostatin significantly inhibited tumorigenesis and 
caused tumor cells apoptosis (Figures 3f and g), which is 
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Figure 1 Myostatin induces metabolic alterations in HeLa cells, (a, b) Effect of myostatin on glucose consumption and lactate production. HeLa cells were treated with 
myostatin (if not specified, at a concentration of 500 ng/ml)^ or PBS for the indicated times, (a) Glucose consumption and (b) Lactate production were monitored, (c) Relative 
lactate level in HeLa, AGZY-83a, MCF-7, SW480 and SH-SY5Y cells cultured in normal DMEM. (d, e) A series of ECAR and OCR in HeLa cells in the presence or absence of 
myostatin were examined by Seahorse XF24 analyzer, (d) ECAR following the addition of glucose defines glycolysis and ECAR following oligomycin represents maximum 
glycolytic capacity. ECAR following treatment with 2-DG represents acidification associated with non-glycolytic activity. Cells are subsequently treated with 4.5 g/l o-glucose, 
1 jM oligomycin and 100 mM 2-deoxyglucose (2-DG). (e) Effect of myostatin on oxygen consumption. OCR prior to the addition of oligomycin (1 fM) and following the 
treatment of FCCP(0.5 /^M) represent the basal mitochondrial respiration and maximal mitochondrial respiration capacity, respectively. OCR following Rotenone (1 ^M) plus 
Antimycin-A(2 /iM) treatment represents non-mitochondrial respiration, (f) Cellular ATP levels in the presence or absence of myostatin. Data are shown as mean 
values ± S.E.M., n = 3. (g) Mitochondrial membrane potential (A^m) was analyzed by flow cytometry with DiOCe^ staining (10 nM). FCCP (2 iM) was used as a positive 
control of A^m loss 



indicated by a much larger total area of cells with condensed 
nuclei (Figure 3h). Together, our results demonstrated that 
myostatin could induce not only cancer cells apoptosis 
in vitro, but also inhibit tumorigenesis in vivo in both 
immune-deficient and immune-competent xenograft tumor 
models. 



Upregulation of VDAC1 is important for myostatin- 
induced apoptosis. VDAC has been recognized as a key 
protein in mitochondria-mediated apoptosis and participates 
in mitochondrial metabolism.^''"^'* Therefore, it is conceivable 
to test the role of VDAC1 in myostatin-induced apoptosis. To 
this end, we first examined several cancer cell lines for 
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Figure 2 Myostatin induces mitochondrial apoptosis in HeLa cells. HeLa cells were treated with myostatin or PBS for the indicated times, (a) Annexin V/PI analysis 
following myostatin treatment. PI, propidium iodide, (b) Time course and dose response of HeLa cells in response to myostatin treatment. The percentage of apoptotic cells 
was detected with annexin V and PI double staining. Data are expressed as mean values ± S.E.M., n = 3. (c) The expression of caspase-9, caspase-3 and PARP were 
examined. Pro-C9, procaspase-9; Pro-C3, procaspase-3, *: non-specific signal, (d) Pan-caspase activity was analyzed using flow cytometry, (e) Clonogenicity of HeLa cells in 
the presence of myostatin. (f, g) Subcellular localization of cytochrome C was examined by western blotting (f) and immunostaining (g) 



VDAC1 expression and sensitivity to myostatin. VDAC1 
expression was markedly induced by myostatin in HeLa, 
AGZY-83a and MCF-7 cells, but remained unchanged in 
SW480 and SH-SY5Y cells (Figure 4a). Interestingly, we 
found that the AGZY-83a and MCF-7 cells, in which VDAC1 
was upregulated by myostatin were also sensitive to 
myostatin-induced apoptosis, whereas the cancer cells in 
which VDAC1 expression was not induced either responded 
to only a limited extent (SW480), or showed complete 
resistance (SH-SY5Y) to myostatin treatment (Figure 4b). 
Notably, VDAC1 induction level and the sensitivity of different 
cancer cell lines to myostatin is positively correlated with their 
glycolytic activities (Figure 1c, measured as lactate produc- 
tion). These data indicate that VDAC1 induction is highly 
correlated with apoptosis induced by myostatin in those 
cancer cell lines addicted to glycolysis, suggesting a 
connection among glucose metabolism, VDAC1 induction 
and cancer cell apoptosis in response to myostatin. 



To determine if the induction of VDAC1 is of functional 
importance for myostatin-induced apoptosis, VDAC1 was 
stably knocked down in HeLa cells. As shown in Figure 4c, 
myostatin-induced apoptosis was substantially attenuated in 
VDAC1 -knockdown strains, m1 and m2. FITC-VAD-FMK 
staining showed significant impairment of myostatin-induced 
caspase activation in cells transfected with VDAC1-shRNA 
(Figure 4d). Together, our data strongly indicate that VDAC1 
plays a critical role in mediating myostatin-induced apoptosis 
of cancer cells. 

The first rate-limiting step in glucose metabolism is the 
conversion of glucose to glucose 6-phosphate (G-6-P) 
catalyzed by hexokinase, which is a key mediator of aerobic 
glycolysis. Hexokinase is known to directly couple intrami- 
tochondrial ATP synthesis to glucose metabolism and to bind 
to VDAC1 and inhibit apoptosis.^^'^^ We then examined the 
expression of Hexokinase II (HKII) and revealed that HKII 
protein level was decreased in HeLa cells treated with 
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Figure 3 Myostatin inhibits tumor growth in xenograft tumor model, (a-e) HeLa cells were injected into nude mice and myostatin were injected intratumorally. (a) Five 
representative xenograft tumors in each group with the treatment of vehicle or myostatin are shown, (b) Tumors weight was measured and shown in the scatter plot (n = 10). 
P= 0.0073. (c) Tumor sizes were measured, and tumor volumes calculated (n = 10 per group), (d) Representative images of the fluorescent TUNEL staining of xenograft 
tumor sections. Arrowheads indicate apoptotic cells. Green: TUNEL-positive. (e) The percentages of apoptotic cells was quantified, and results are presented as mean 
values ±S.E.M. of 10 sections. At least 200 cells from each section were counted, (f-h) B16F10 cells were injected into C57BL/6 mice and myostatin were injected 
intratumorally every other day. (f) Seven representative xenograft tumors in each group with the treatment of vehicle or myostatin are shown, (g) Tumors weight was measured 
and shown in the scatter plot 13). P= 0.0055. (h) Representative images of the hematoxylin and eosin staining of xenograft tumor sections. Arrowheads indicate the 
regions enriched with apoptotic cells, indicated by condensed nuclei 



myostatin in a time-dependent manner (Figure 4e). Fractiona- 
tion analysis further showed that mitochondrial HKII, consis- 
tent with the pattern of total HKII, decreased at 24 h treatment 
(Figure 4f). This is in contrast to the expression level of 
VDAC1 , which was found to be highly induced in the presence 
of myostatin in a time-dependent manner (Figures 4a, e and f). 
However, other mitochondrial markers, including Tim23 and 
UCP2 for inner membrane, and cyclophilin D and MnSOD for 
matrix proteins remained unchanged. To further explore 
whether HKII is involved in protecting cells from myostatin- 
induced apoptosis, HeLa cells with stably overexpressed HKII 
were treated with myostatin, and both Annexin V/PI staining 
and FITC-VAD-FMK showed a significant reduction of 
apoptosis in those cells with overexpressed HKII (Figures 
4g and h). In addition, we were able to show myostatin- 
induced expression of VDAC1 was due to neither 



transcriptional regulation nor the increased de novo protein 
synthesis (Supplementary Figure S6). Taken together, our 
data revealed that the observed changes in HKII and VDAC 
protein levels are myostatin-specific responses, and that 
these changes are closely linked with myostatin-induced 
apoptosis in cancer cells. 

ATP depletion is responsible for myostatin-induced 
VDAC1 upregulation and apoptosis. We next are 
prompted to ask if metabolic changes and ATP depletion is 
causally linked with mitochondrial-dependent apoptosis in 
cancer cells. To test this, we used 2-deoxy-glucose (2DG) to 
mimic ATP depletion (Figure 5a) through glucose depriva- 
tion. Western blotting showed that 2DG alone was sufficient 
to induce VDAC1 expression (Figure 5b). HeLa cells were 
then treated with myostatin in the presence of 2DG and we 
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Figure 4 VDAC1 is important in myostatin-induced apoptosis. (a) Expression of VDAC1 in different cell lines with the treatment of myostatin for the indicated times, 
(b) Different cancer cell lines demonstrate varied sensitivity to myostatin treatment. Apoptosis was analyzed by Annexin V/PI staining, (c, d) HeLa cells stably transfected with 
scrambled or VDAC1 shRNA were treated with myostatin or PBS for 36 h at the indicated concentrations. siVD-ml : stable clone 1 ; siVD-m2: stable clone 2. Apoptosis (c) and 
pan-caspase activity (d) were analyzed using flow cytometry, (e) HeLa cells were treated with myostatin and cell lysates were analyzed by western blotting with antibodies 
for VDAC1, HKII, Cyclophilin D, UCP2, Tim and MnSOD. (f) Localization of HKII, VDAC1 and Tim in different fractionations of HeLa cells with the treatment of myostatin. 
(g, h) HeLa cells, stably transfected with myc-HKII (upper panel of g), were treated with PBS or myostatin for the time as indicated. Apoptosis (lower panel of g) and pan- 
caspase activity (h) were analyzed using flow cytometry 



found that 2DG significantly enlianced myostatin-induced 
apoptosis of the HeLa cells (Figure 5c). Mitochondrial 
respiratory inhibitors such as rotenone (complex I), antimycin 
A (complex III) and oligomycin (complex V) inhibited 
OXPHOS resulting in ATP depletion, and interestingly, they 
all increased VDAC protein expression (Supplementary 



Figure S7). The effects of mitochondrial inhibitors further 
support the observation that VDAC1 can be induced in 
response to ATP depletion and consequentially lead to 
apoptosis. Interestingly, inosine (IN) and sodium pyruvate 
(SP), which are routinely used to restore cellular ATP 
levels, ^^'^^ rescue ATP depletion and partially inhibited 
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Figure 5 ATP depletion mediates myostatin-induced VDAC1 upregulation, which contributes to apoptosis. (a, b) HeLa cells cultured in low glucose (5.5 mM) DMEM were 
treated with vehicle, 8 mM, 16 mM or 32 mM 2DG for 24 h. Cells were harvested for the ATP assay (a) and western blotting analysis (b). (c) HeLa cells pre-incubated with or 
without 2DG (32 mM) for 30 min were treated with PBS or myostatin for 24 h, followed by Annexin V/PI staining, (d) HeLa cells were incubated in glucose-free DMEM in the 
presence or absence of sodium pyruvate (1 0 mM) or inosine (20 mM), followed by treatment with myostatin (500 ng/ml) for 6 h. Cellular ATP levels were examined (n = 3). 
(e) HeLa cells were treated as in (d), and cells lysates were subjected to western blotting analysis with anti-VDAC1 antibody, (f) HeLa cells were treated as in (d) with extended 
myostatin treatment for 12 h, followed by Annexin V and PI double staining 



myostatin-induced VDAC1 upregulation (Figure 5e) and 
apoptosis (Figure 5f). Combined, all of these results support 
the conclusion that ATP depletion is responsible for 
myostatin-induced VDAC1 upregulation and apoptosis. 

Mitochondrial translocation of Bax is involved in 
myostatin-induced apoptosis. Previous studies have 
shown that detachment of HK II from the mitochondria 
promoted the translocation of Bax to mitochondria for 
subsequent cytochrome c release.^^'^^'^° Consistent with 
this, we found that myostatin treatment causes Bax 
translocation to mitochondria as revealed by immunostaining 
of Bax protein (Figure 6a). This observation was further 
confirmed by fractionation analysis as shown in Figure 6b. In 
addition, the expression profile of Bcl-2 family members 
following myostatin treatment clearly showed that the anti- 
apoptotic Bcl-2 family proteins Bcl-2, Mcl-1 remain unaltered 
(Figure 6c). The protein levels of the pro-apoptotic proteins of 
Bax and Bak are also not changed, while Bim and Bik levels 
were reduced following myostatin treatment (Figure 6c). 
Further analysis shows that myostatin-induced apoptosis is 
largely inhibited in the Bax-deficient HCT116 cells 
(Supplementary Figure S8A), indicating the specific role of 
Bax in myostatin-induced apoptosis. Although we did not 
detect direct interaction between Bax and VDAC1 (likely 
because Bax:VDAC1 interaction is transient and weak), we 
observed significant inhibition of Bax translocation and 
cytochrome c release when VDAC1 was silenced with 
shRNA (Figures 6d and e). As expected, overexpression of 



Bcl-2 also prevents myostatin-induced Bax translocation and 
apoptosis (Supplementary Figures S8B-D). We thus con- 
clude that chronic exposure of myostatin causes metabolic 
stresses leading to VDAC1 upregulation, Bax translocation to 
mitochondria and subsequent apoptosis in cancer cells. 

Discussion 

Within the present study we have provided compelling 
evidence that myostatin shifts glucose metabolism towards 
glycolysis and induces apoptosis in cancer cells. Exposure to 
myostatin resulted in mitochondrial morphological, physiolo- 
gical and functional changes, leading to the mitochondria- 
dependent apoptosis in cancer cells (Figure 6f). We found that 
myostatin affects the protein levels of HK II and VDAC1, two 
key regulators of glucose metabolisms. The alteration of 
mitochondrial metabolism could lead to dissociation of HKII 
from the mitochondria and Bax translocation to the mitochon- 
dria, resulting in the metabolic stress-induced apoptosis. 

Although the importance of activin/TGF-j5 and their recep- 
tors on apoptosis and cancer biology is well recognized, to our 
best knowledge, our work represents the first report showing 
that myostatin induces mitochondria-dependent apoptosis in 
cancer cells. TGF-p and related family members may 
transcriptionally regulate apoptosis-related genes, such as 
XIAP, Bcl-2/Bcl-xl or Bim. Alternatively, they regulate apop- 
tosis through non-canonical pathways, such as JNK and Akt 
signaling, leading to the phosphorylation and activation of the 
pro-apoptotic molecules.^'' '^^ Nonetheless, upon myostatin 
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Figure 6 Mitochondrial translocation of Bax is involved in VDAC1 -mediated apoptosis, which is induced by myostatin. (a) Left panel: HeLa cells treated with myostatin for 
36 h were fixed and stained with anti-BAX (green) and anti-TIIVI (red) antibodies. Right panel: Statistics of HeLa cells with redistributed Bax in the presence of myostatin for 
36 h. The data were shown as mean ± S.E.IVI. (b) Subcellular localization of VDAC1 and Bax were examined by western blotting. Tim was shown as a loading control of 
mitochondrial sample; actin was shown as an loading control of cytosolic fraction and whole cell lysate. (c) Expression profile of major members of Bcl-2 family, (d) Scramble 
and siVD-ml cells were treated with PBS or myostatin for 36 h. Immunofluorescence was performed with antibodies for Bax (green) and Tim (red), (e) HeLa cells were treated 
as in (d). Green: Cytochrome C; Red: Tim-Cy3. (f) A model for the mechanism underlying myostatin induction of VDAC1 -mediated apoptosis in cancer cells 



treatment, we did not find any induction or phosphorylation of 
pro-apoptotic Bcl-2 family proteins, such as Bim and Bad, nor 
did we see the drastic changes of anti-apoptotic Bcl-2 
molecules. Although Bcl-xl was slightly decreased in 
response to myostatin, the functional analysis did not show 



any correlation with myostatin-induced apoptosis. Instead, we 
clearly showed that myostatin activates Bax-dependent 
apoptosis. Bax translocation to the mitochondria upon 
myostatin treatment is highly correlated with the onset of 
apoptosis and cytochrome c release. 
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It is recognized that dissociation of hexokinase from 
mitochondrial VDAC regulates Bax translocation to the 
mitochondria.^^'^^'^°'^^ Our finding is the first to suggest that 
such a mechanism is likely involved in myostatin-induced 
apoptosis. We demonstrated that VDAC1 plays a critical role 
in myostatin-induced apoptosis of cancer cells (Figures 4 
and 5). One theory, albeit controversial,^"^ is that closure of the 
VDAC1 channel prevents efficient exchange of ATP and ADP 
between the cytosol and mitochondrial matrix, leading to 
mitochondria swelling and outer membrane permeabilization 
(MOMP), and thereby release of pro-apoptotic factors from 
the intermembrane space. ^^'^^'^^'^^ Indeed, we observed a 
loss of mitochondrial membrane potential, indicative of the 
opening of mitochondrial permeability transition pore (mPTP), 
during the onset of myostatin-induced apoptosis. Interest- 
ingly, we also found that Gyclosporin-A (CsA), a selective 
inhibitor of mPTP (Supplementary Figure S9), can compro- 
mise the apoptosis-inducing function of myostatin, implying 
that mPTP may also contribute to myostatin-induced apopto- 
sis. Alternatively, VDAC could interact with HKII or Bcl-2 
family proteins, thus regulating apoptosis.^^'^^ Indeed, we 
observed that dissociation of HKII from the mitochondria 
occurs, concomitantly, with the translocation of Bax to the 
mitochondria. Akt and GSK3j5 are implicated in regulating the 
association of HKII with VDAG1 Interestingly, we 

previously reported that myostatin could activate GSK3j5 
and induce cell-cycle arrest through the inhibition of Akt in the 
myoblast,^ which aids in understanding the case in cancer 
cells, in which prolonged treatment of myostatin may not 
sustain Akt activity, and instead leads to the dissociation of 
HKII from mitochondria. 

It is intriguing that the changes in VDAG and HKII protein 
levels are inversely correlated upon myostatin treatment. As 
was reported, depending on the cellular metabolic state, HKII 
may either promote or inhibit apoptosis.'*^ Interestingly, we 
found that perturbation of intracellular ATP levels, either by 
application of 2DG or inhibitors of respiratory chain, could 
induce VDAG1 protein levels as we observed in myostatin- 
treated cancer cell, whereas replenishment of ATP prevented 
VDAG1 upregulation and myostatin-induced apoptosis. It is 
possible that, when ATP level are decreased, as a compen- 
satory response the cell may enhance the stability of VDAG1 
for better transportation of ATP to supply enough energy for 
cell survival. Apart from GSK3j5, it is known that a number of 
metabolism-related kinases, such as PKA and PKG phos- 
phorylate VDAG1 and regulate the interaction between 
VDAG1 and Bcl-2 family proteins, which in turn regulates 
respiration and apoptosis. ^^'^^ Additionally, it was known that 
phosphorylation of VDAG1 may trigger conformational 
changes of VDAG1 and inhibit its degradation, resulting in 
VDAG1 accumulation.^^'"^^ Indeed, the real-time PGR results 
showed that myostatin did not alter the mRNA level of VDAG1 , 
whereas induction of VDAG1 in response to myostatin 
treatment could not be blocked by inhibiting de novo protein 
synthesis with cycloheximide (Supplementary Figure S6D), 
which suggests that myostatin could possibly enhance 
VDAG1 stability. On the other hand, phosphorylation of 
VDAG1 resulted in the dissociation of HK from the mitochon- 
dria,^^ thus reducing its stability in cytosol. Intriguingly, the 
total HKII activity maintained to the same level with untreated 



cells, even though we observed drastic decline of its 
expression at 24 h treatment (Supplementary Figure S10). 
While the exact mechanism of VDAG1 and HKII stability 
awaits further investigation, our data suggest that altered ATP 
levels and metabolic stress affect mitochondrial-dependent 
apoptosis. 

Myostatin is regarded as a potential molecular target for a 
number of diseases, including muscular dystrophy, sarcope- 
nia, cachexia and diabetes.'^^'^'^ Our results are the first to 
suggest that myostatin could be useful in fighting cancers that 
have altered glycolytic activity and VDAG1 response. Our 
results from both in vitro and xenograft tumor models further 
indicate that focal expression or intratumoral injection of 
myostatin may be effective ways to combat cancer, while 
avoiding deteriorating effects on muscle. Both HKII and 
VDAG1 could serve as targets for cancer therapy as 
previously suggested. ^^ "^^""^^ As a matter of fact, myostatin 
is known to be secreted in a latent form and is maintained at a 
comparatively high concentration in muscle, where it tends to 
be activated under acidic and hypoxic conditions due to high 
levels of glycolysis. ^^'"^^ Therefore, it is not surprising that 
myostatin might have specific proapoptotic functions in tumor 
tissue, but not in normal tissue, because most tumor cells 
predominantly produce ATP through glycolysis, and thus form 
a heterogeneous acidic and hypoxic microenvironment, which 
is similar to the case of muscle after exercise. In order to testify 
this hypothesis, we designed a series examination on different 
steps of myostatin-induced apoptosis in different cancer cell 
lines. Surprisingly, the cell-type-dependent sensitivity to 
myostatin significantly correlates to their glycolytic activity 
(r= 0.991 , P= 0.0005) (Supplementary Figures S11 A and B; 
Figure 1c), whereas shows no correlation with the basal 
expression level of AGTRIIB, VDAG1 (Supplementary 
Figures S11G-E), or the quiescent ATP level 
(Supplementary Figures S11H and I) in different cells, 
although the decrease of ATP and the induction of VDAG1 
have been shown to be necessary for myostatin-induced 
apoptosis. Moreover, Lactate accumulation in the medium 
could be influenced by several factors besides the real 
glycolysis level, such as the proliferation rates of different 
cancer cells. However, we do see the big variation among 
cells derived from different cancers, but they showed no 
correlation with myostatin susceptibility (Supplementary 
Figures SI IF and G). Together, our results may help to 
understand, at least in part, why there is a low incidence of 
cancer in skeletal muscle, and how physical exercise could 
prevent cancer. During exercise, glucose is rapidly utilized by 
muscle and lactate is highly produced, as observed in our 
system. These findings may also aid in understanding why 
cells addicted to glucose are more sensitive to myostatin, and 
serve to explain why certain malignant tumor cells are more 
sensitive to myostatin than cells in normal tissue. 

Materials and Methods 

Cell culture and transfection. HeLa, AGZY-83a, MCF-7, SW480, SH- 
SY5Y and B16F10 cells were grown in DMEM supplemented with 10% FBS 
(Hyclone, Logan, UT, USA) and 1% penicillin-streptomycin at 37 °C under 5% 
CO2. The targeting sequence in VDAC1 for RNA interference, 5' aaAGT- 
GACGGGCAGTCTGGAA 3', was reported earlier by Yuan et al."^^ Primers for 
shRNA were designed according to the manufacturer's instructions (Ambion, 
Carlsbad, CA, USA), and cloned into psilencer2.1-neo vector. To establish stable 
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cell lines, HeLa cells were transfected with sh-VDAC1 or pcDNA4T0B-myc-HKII, 
and screened with neomycin and ziocin, respectively. Recombinant myostatin was 
expressed and purified as described previously.^^ Exponentially growing cells 
were subjected to various treatments. 

Reagents and antibodies. Antibodies against PARP, cytochrome c and 
Tim were obtained from BD Biosciences (San Jose, CA, USA). Caspase-9 and 
caspase-3 antibodies were acquired from Cell Signaling (Danvers, MA, USA), and 
the anti-VDAC1 monoclonal antibody from Abeam (Cambridge, MA, USA). SP, 
inosine and 2-deoxy-D-glucose were purchased from Sigma Chemical Co. 
(St. Louis, MO, USA). TMRM and DiOCe^ were acquired from Invitrogen 
(Carlsbad, CA, USA). The caspACE FITC-VAD-FMK in situ marker was purchased 
from Promega (Madison, Wl, USA). 

SDS-PAGE and western blot. Cells or membrane fractions were lysed in 
lysis buffer (10 mM HEPES, pH 7.4, 2mM EGTA, 0.5% NP-40, protease 
inhibitors). Equivalent protein quantities (20 ^g) were subjected to SDS-PAGE, 
and transferred to nitrocellulose membranes. Membranes were probed with the 
indicated primary antibodies, followed by the appropriate HRP-conjugated 
secondary antibodies (KPL, Gaithersburg, MD, USA). Immunoreactive bands 
were visualized with a chemiluminescence kit (Pierce, Rockford, IL, USA). 

Subcellular fractionation. Subcellular fractionation was performed as 
previously described.^" Cells were harvested and resuspended in hypotonic buffer. 
After gentle homogenization with a Dounce homogenizer, cell lysates were 
subjected to differential centrifugation. The resulting membrane fractions were 
lysed, and analyzed using western blot. 

Immunofluorescence microscopy. Cells were grown to 60% confluence 
on a coverslip. After treatment, cells were washed twice with PBS, and fixed with 
freshly prepared 3.7% formaldehyde at 37 °C for 15 min. Antigen accessibility was 
increased by treatment with 0.2% Triton X-100. Cells were incubated with primary 
antibody for 1 h, and, after washing, stained with a secondary antibody for a 
further 45 min. Cell images were captured with an LSM 510 Zeiss confocal 
microscope (Carl Zeiss, Jena, Germany). DAP! staining was used to identify 
nuclei. 

Apoptosis assay. Apoptosis was assessed by detecting phosphatidylserine 
exposure on the cell membrane with Annexin V staining as described previously.^'' 
In situ caspase activation was determined by detecting FITC-VAD-FMK 
fluorescence with flow cytometry as described previously.^^ TUNEL analysis 
was performed with an in situ death detection kit (Roche Molecular Biochemicals, 
Mannheim, Germany) as described previously.^^ 

Clonogenic assay. In brief, cells were pre-incubated for various periods of 
time in the presence of myostatin. Subsequently, cells with each treatment were 
plated in a six-well plate at 1000 cells per well in triplicates. Cell colonies were 
stained with crystal violet and then counted 10 days after initial plating. 

Mitochondrial membrane potential assay and metabolic events 
measurements. Mitochondrial membrane potential was measured with 
DiOCe^ as described previously.^^ Metabolic events were detected according to 
the protocols from Methods in Cell Biology.^^ Intracellular ATP levels were 
determined using a commercially available colorimetric ATP assay kit (Biyuntian 
Co., Shanghai, China) according to the manufacturer's instructions. ATP levels 
were normalized to protein levels. Glucose consumption and the lactate level were 
evaluated with commercial kits (Biovision, Mountain View, CA, USA). 

OCR and ECAR were evaluated by using a Seahorse XF24 Analyzer (Seahorse 
Bioscience, North Billerica, MA, USA) as described.^"^ In brief, 2x10"^ HeLa cells 
were seeded 24 h before myostatin treatment. The OCR is measured sequentially 
before or after the addition of oligomycin (1 fiM), FCCP (0.5 ^M) and rotenone 
(1 i^iU) plus antimycin-A (2 ^M). The ECAR is determined sequentially before or 
after the addition of o-glucose (4.5 g/l), oligomycin (1 ^M) and 2-DG (100 mM). 

Tumor xenograft studies. Five-week-old female nu/nu nude mice were 
injected in the right posterior flank region with 2.5 x 10^ HeLa cells in 0.1 ml of 
medium. Seven days after implantation, when most of the tumors had grown up to 
75 mm^, animals were randomly selected to receive either vehicle (n = 10) or 
treatment (n= 10), followed by direct intratumoral injection of PBS or myostatin at 



a concentration of 0.5 mg/kg/day. Tumor sizes and body weights of mice were 
recorded every other day. Tumor dimensions were measured with calipers, and 
the volumes calculated using the following formula: volume = 0.5 x / x w^, with / 
representing the maximal length and w the width. After treatment for 3 weeks, all 
mice were sacrificed and carcinoma tissue evaluated using fluorescence- 
conjugated TUNEL assay. 

For melanoma xenograft model, 1x10^ B16F10 cells were injected subcuta- 
neously into C57BL/6 mice. Ten days after implantation, when palpable tumor formed, 
myostatin were injected into the tumor mass at a dose of 0.5 mg/kg every other day. 
PBS-injected mice served as controls (n = 13). Tumors were measured every 2 days, 
and volume was calculated (length x width^ x 0.5). All mice were sacrificed 21 days 
after cell injection and tumor tissue were collected and weighed. A 5-/im thick formalin- 
fixed tumor sections were stained with hematoxylin and eosin. 

Statistical analysis. In quantitative analyses using cultured cells represented 
as histograms, values were obtained from three independent experiments, and 
expressed as the mean±S.E.M. Statistical analysis was performed using the 
Student's f-test, with P-values<0.05 considered significant. *P<0.05 and 
**P<0.01 versus the corresponding controls are indicated. All statistical data 
were calculated with GraphPad Prism software (GraphPad Software, Inc., La Jolla, 
CA, USA). 
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